Liquid transport generated by a flashing field-induced wettability ratchet 
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Abstract 

We develop and analyze a model for ratchet-driven macroscopic transport of a continuous phase. The 
transport relies on a field-induced dewetting-spreading cycle of a liquid film with a free surface based on 
a switchable, spatially asymmetric, periodic interaction of the liquid-gas interface and the substrate. The 
concept is exemplified using an evolution equation for a dielectric liquid film under an inhomogeneous 
voltage. We analyse the influence of the various phases of the ratchet cycle on the transport properties. 
Conditions for maximal transport and the efficiency of transport under load ai^e discussed. 
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Brownian ratchets present a well established way to induce directed motion of particles in spa- 
tially extended systems without global gradients, i.e. without a globally broken spatial symmetry 
] examples include colloidal particles suspended in solution that move when exposed 
to a sawtooth dielectric potential [41] and the selective particle filter formed by a micro-fabricated 
silicon membrane with asymmetrical bottleneck- like pores under application of an oscillating pres- 
sure gradient [jsO . Brownian ratchets are also thought to represent the underlying mechanism of 
molecular motors that are responsible for the active transport of molecules along filaments in bi- 
ological cells [6y. Brownian transport is based on a principle first pointed out by Pierre Curie [|70 
stating that also a macroscopically symmetric constellation may induce macroscopic transport if it 
exhibits local asymmetries, e.g. a periodic asymmetric potential that varies on a small length scale. 
However, the system has to be kept out of equilibrium, for instance, by a chemical reaction an 
oscillating pressure [|5D or electric potential j^. 

Ratchets are not only used to transport or filter discrete objects like colloidal particles or macro- 
molecules. A ratchet geometry may also serve to induce a macroscopic transport of a continuous 
phase using local gradients only. Experimentally investigated examples include a secondary liq- 
uid flow triggered in Marangoni-Benard convection involving a solid substrate with asymmetric 



lOb. Micro- 



grooves [Sj. On a similar substrate Leidenfrost drops perform a directed motion |19|, 
drops confined in asymmetrically structured geometries move when vibrating the substrate or ap- 
plying an on/off electric field [lllll . 

In this Letter we propose and analyse a simple generic model for ratchet driven free-surface 
flows resulting in the macroscopic transport of a continuous phase. The driving flashing ratchet is 
based on a switchable spatially periodic but asymmetric interaction of the liquid-gas interface and 
the solid substrate like, for instance, a field-induced wettability. One can envision several ways 
to experimentally realise the spatially inhomogeneous interaction and the switching in time. A 
simple setup consists of a thin film of dielectric liquid in a capacitor producing a sawtooth voltage 
profile when switched on. The electro-dewetting provoked by the application of a homogeneous 



electrical field is already used to structure thin films of liquid polymers [|12LI13L 

An idealized electrical wettability ratchet (sketched in Fig. [U works as follows. We consider 
a flat wetting film stable in the absence of an electric field (Fig. [I](a)). Upon switching on an 
electric field at t = (Fig. [11(e)) the film is destabilized by both, the overall electric field and 
ocal gradients (Fig.[T](d)). This behavior is similar to dewetting on a heterogeneous substrate 



Its 



160. After a transient it results in the collection of all the liquid in drops situated at the patches 
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FIG. 1 : Panels (a) to (c) illustrate the working principle of a fluidic ratchet based on a switchable wettability 
that causes dewetting-spreading cycles, (d) illustrates the spatial asymmetric periodic voltage profile <l>(x) 
responsible for the wettability pattern and (e) indicates the time dependence Q(t) of the switching in relation 
to the dewetting and spreading phases in (a) to (c). 

of maximal voltage (Fig.[T](b)). After switching off the field at t = r (Fig.[T](e)) the drops spread 
into a homogeneous wetting layer (Fig. [I](c)). Finally, at t = T the field is switched on again and 
a new cycle starts. A typical realistic film dynamics during one ratchet cycle is shown in Fig. [2l 
For the geometry in Fig. [11(d), with a voltage profile shaped like a sawtooth skewed to the right, 
each cycle transports liquid to the right. 

The inhomogeneous electric field represents one out of several possible choices for a switchable 
spatially asymmetric interaction of liquid-gas interface and substrate. Others not pursued here 
include switchable brushes lIlTn . heating |ll8i[l9(] or the switching of optical substrate properties, 
i.e. van der Waals interactions. Note however, that although length and time scales and the specific 
form of the resulting interaction potentials will differ, in principle all these realizations can be 
mapped onto our model using appropriate pressure terms corresponding to an effective switchable 
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'wettability'. 

In presenting the dynamical model we restrict our attention to a two-dimensional system corre- 
sponding to a shallow open channel geometry and neglect the influence of the channel walls. The 
interplay of field-driven dewetting, liquid motion caused by local lateral gradients and spreading 
for small scale systems, i.e. of wettability and capillarity dominated fluid flow, is well described 
using an evolution equation for the film thickness profile fllOfl. A derivation from the basic trans- 
port equations using lubrication approximation gives 



dfh = 



3r] 



(d^p- /cxt) 



(1) 



where 7 and 77 are the surface tension and dynamic viscosity of the liquid, respectively, /cxt 
denotes an external force in the direction of the x-axis. The change in time of the film thickness 
profile is given by the divergence of the flow, expressed as the product of a mobility and a pressure 
gradient and/or external force. The velocity field within the film is fully determined by the film 
thickness profile. The velocity parallel to the substrate is u{x,z) = — zh) d^p and the 

vertical component is obtained using continuity. The pressure 



P 



■'jd^^h - n(/i,x,t) 



(2) 



contains the curvature pressure —'jdxxh and the disjoining pressure Il{h,x,t). The latter com- 
prises the effective interactions between the liquid-gas interface and the substrate, i.e. the wetta- 



bility properties 



21, 



2211 . Note, that the lubrication approximation can be formally applied 



to systems involving small surface slopes only. However, even for partially wetting systems with 
large contact angles, the lubrication approximation predicts in most cases the correct qualitative 
behavior [ 2o|l . 

As model system we use a dielectric oil in a capacitor of gap width d. The oil wets the lower 
plate and does not wet the upper plate corresponding to the van der Waals disjoining pressure 
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with the Hamaker constants Ai > and A„ < 0. A dielectric fi 



voltage Uo is subject to an electrical 'disjoining' pressure u2, 



(3) 



m in a capacitor with an applied 



(4) 
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where Eq and si are the absolute and relative dielectric constant, respectively and e = Exje^ 
denotes the ratio of the relative dielectric constants of the liquid and the gas phase. The modulation 



of the electric field is periodic in x, but breaks parity, i.e. there exists no reflection symmetry in x. 
This ratchet potential is periodically switched on and off. In consequence 11 in ([21) takes the form 

n(/i, X, t) = n{t) Uei{h) + U,dw{h). (5) 

where the temporal variation and spatial variation are defined in Fig. [T](a) and[T](b), 
respectively, with (l/L) dx = 1. 




FIG. 2: Space-time plot of one temporal period of the evolution of the film thickness profile for a film of 
dielectric liquid in a capacitor. Shown is one spatial period. The different phases (i)-(iv) indicated to the 
right are explained in the main text. Parameters aieh = 0.5, tp = 0.5, L = 32, = 5, T = 5000, w = 1, 
A = 0.001, e = 2.5 and /ext = 0. The starting time is well after initial transients have decayed. 

To obtain a minimal set of parameters we introduce the scales S^rj/dKli, ^y^yd/K^i, and d for 
t, X, and h, respectively. The electrostatic 'spreading coefficient' is defined by Kei = eoEiUQ /2d'^. 
The resulting dimensionless equations correspond to ([U)-© with St] = 'j = d = eo^i^^o/^ = 
1. For simplicity we assume Au = —Ai and obtain the dimensionless Hamaker constant A = 
Ai/Qud^Kei- All results are given in terms of dimensionless quantities. 

We characterize the ratchet using two measures, the asymmetry ratio = (L — a) /a of the 
spatial variation $(x) and the flashing ratio = r/(T — r) of the temporal switching Vt{t). Zero 
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net transport is expected for a symmetrical ratchet, i.e. for = 1. The resulting transport is 
quantified by the mean velocity v = j/h where j = (l/TL) dt dxj{x, t) is the mean flow 
along the substrate and t) = —h^{dxP — /cxt)- 

Figure |2] shows a typical example of the film evolution during one ratchet cycle allowing to dis- 
tinguish four phases, (i) When the ratchet is switched on the film is nearly flat but rapidly evolves 
a surface instability with a wavelength given approximately by the corresponding spinodal length 
boll , (ii) Next, the evolving profile coarsens accelerated by the gradients of the ratchet potential. 



(iii) In an ideal situation only one drop remains, corresponding to the equilibrium structure on the 
heterogeneous wettability pattern produced by the ratchet potential (iv) Finally, after switch- 
ing off the ratchet the drop starts to spread rapidly under the influence of van der Waals forces 
until the next cycle starts. 




FIG. 3: Dependence of the mean velocity v on the liquid level in the capacitor, i.e. on the mean film 
thickness h, for various flashing periods T as indicated in the legend. The remaining parameters are as in 
Fig.lH 



The competing influence of the various parameters allows to tune the transport properties by 
adjusting the relative importance of the various phases. As shown in Fig.[3]the transport is strongest 
for intermediate film thicknesses. A film of intermediate thickness in a homogeneous electric field 
= = 1) dewets spinodally with a wavelength well below the used spatial period of the 



ratchet. Therefore the ongoing coarsening interacts with the flow induced by the spatial asymmetry 
and leads to a nontrivial dependence of the mean velocity on the film thickness. For very small or 
large film thicknesses the stabilizing van der Waals terms are dominating and the mean velocity 
approaches rapidly zero. Note also the flow reversal, albeit with a small mean velocity, for large 
h > 0.85. Further calculations (not shown) demonstrate the monotonous increase of the mean 
velocity with increasing asymmetry ratio or amplitude ip of the ratchet potential for otherwise 
constant parameters. 
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FIG. 4: Variation of the mean velocity v in dependence of (a) the flashing ratio lo (for T=500) and (b) the 
flashing period T (for co = 1) for various film thicknesses as given in the legends. The remaining parameters 
are as in Fig. |2] 

The flashing characteristics of the ratchet have a very pronounced effect on the transport. Fig- 
ures |4](a) and (b) present the non-monotonous dependencies of the flow on the flashing ratio to and 
period T, respectively. For small u the flow is practically zero, since the time for dewetting is too 
short to trap a considerable amount of liquid at the spots of high wettability. Increasing u increases 
the velocity until it reaches a maximum at ^ 1. Beyond the maximum the flow decreases again 
because less and less time remains for the spreading. For large u one observes a flow reversal. 

The dependence of the velocity on the flashing period T is similar but shows around the flow 
maximum a particularly interesting non-monotonous behavior that is due to coarsening. For small 
periods the fluid has neither enough time to dewet nor to spread and the resulting mean transport 
is small. For large periods both processes reach the respective equilibrium structure well before 
the next switching, i.e. most time is spend waiting and the mean velocity decreases approximately 
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as I IT. 

Assuming an ideal combination of ratchet properties and intrinsic length and time scales, and 
approximating the intermediate drops at the field maxima as point-like objects, all the liquid moves 
by a — L/2 in one cycle of length T. For a film of mean thickness li this yields a mean flow of 
jideai = h{a — L/2)/T, i.e. the mean velocity is i)ideai = {a — L/2)/T. The maximal spatial 
asymmetry is given by a = L, i.e. jmax = hL/2T and Vmax = L/2T. In a real system, however, 
several effects keep the mean velocity below the ideal value. For optimal parameters in simulations 
we reach about 30% of Videai (Fig. |4](b) for h = 0.7 and T ^ 10000). 

For practical applications transport against an external force is of crucial importance. In the 
reference system of a uniform substrate no work is performed and all energy is lost via viscous 
dissipation. This is not the case for a thin film ratchet under load, i.e. /cxt 7^ in ([T])- Switching 
off the ratchet, /cxt < induces a negative macroscopic flow. Simple mechanical realizations of 
/cxt are an inclined substrate or centrifugal forces yielding a constant /cxt independent of the film 
thickness. The flashing ratchet generates a positive flow, i.e. mechanical work is performed against 
the external force. The transport under load is characterized by the (mechanical) energy transport 
efficiency 

W 

where ^ ^ 

W = ^ [ dx [ dtjUt (7) 



denotes the mechanical work performed per unit time and 



denotes the interaction energy consumed per time. The mean velocity and transport efficiency are 
shown in Figs.[5](a) and (b). The transport efficiency behaves non-monotonically. It increases for 
small loads, reaches a maximum and decreases for higher loads. It becomes negative when the 
flow reverses. Then the ratchet can no longer perform work against the external force. 

Finally, we discuss the feasibility of a fluid ratchet employing spatially inhomogeneous electric 
fields by estimating the acting forces and relevant time and length scales for millimetric, micro- 
scopic and nanoscopic fluidic systems. Typical material constants for thin films in capacitors are 



taken from ||23, 



25|1 and references therein. Applying a voltage of lOOV over a gap of 2 mm 
width, a 1 mm film of silicon oil feels an electrostatic pressure of Pei ~ O.lNm"^ corresponding 
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FIG. 5: The influence of an external load /cxt- Shown are the variations of (a) the mean velocity v and (b) 
the transport efficiency VcS with /cxt for various mean film thicknesses h as indicated in the legends. The 
remaining parameters are as in Fig.|2]except for T = 500. 

approximately to the curvature pressure in a drop of 1 mm height with = 10 ° equilibrium con- 
tact angle. Thinner films can be influenced by a lower voltage lll2ll . For instance, for a 10 fim oil 
film in a 20 /^m gap for f/g = 10 V we find Pei ~ 10Nm~^, equivalent to the curvature pressure 
in a droplet of 10 fxm height and 9q = 10°. In principle, nanofluidic transport of polymer melts 
is also feasible. For a 30 nm film of liquid polystyrene in a 100 nm gap which is electrostatically 



patterned (Uq = 10 V) [|25|] . pei ~ lO'^N/m^ equivalent to the curvature pressure for a droplet of 
30 nm height and 300 nm width. Note, however, that for ultrathin films the disjoining pressure is 
of the same order of magnitude [|22n . 

Beside the length scales also the time scales are of uttermost importance. The scale r = 
3'yr]/dKli reflects the relevant properties responsible for the relaxation towards the flat film. For 
a 5cS silicon oil the estimate gives a timescale of 10 to 10^ s, i.e. the viscous flow is rather slow. 
Using, however, water films of thicknesses between lO/im (Uq = lOV) and 1 mm (Uq = lOOV) 
the time scale ranges from 10^^ to 10"'^ s resulting in fast transport. 

In conclusion, we have shown that a flashing ratchet produces a macroscopic transport in a 
liquid film with a free surface on a feasible time scale. There exist regimes of maximum transport 
selected by the spatial and temporal properties of the ratchet, which depend on the characteristics 
of the thin film. For a spatial period of the ratchet which is considerably larger than the length scale 
of the spinodal instability of the flat film subjected to a homogeneous potential the coarsening 
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dynamics influences the transport in a non-trivial way. 
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